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ABSTRACT 

The secular evolution of the purely general relativistic low angular momentum accretion flow 
around a spinning black hole is shown to exhibit hysteresis effects. This confirms that a sta- 
tionary shock is an integral part of such an accretion disc in the Kerr metric. The equations 
describing the space gradient of the dynamical flow velocity of the accreting matter have been 
shown to be equivalent to a first order autonomous dynamical systems. Fixed point analy- 
sis ensures that such flow must be multi-transonic for certain astrophysically relevant initial 
boundary conditions. Contrary to the existing consensus in the literature, the critical points 
and the sonic points are proved not to be isomorphic in general, they can form in a com- 
pletely different length scales. Physically acceptable global transonic solutions must produce 
odd number of critical points. Homoclinic orbits for the flow flow possessing multiple criti- 
cal points select the critical point with the higher entropy accretion rate, confirming that the 
entropy accretion rate is the degeneracy removing agent in the system. However, heteroclinic 
orbits are also observed for some special situation, where both the saddle type critical points 
of the flow configuration possesses identical entropy accretion rate. Topologies with hetero- 
clinic orbits are thus the only allowed non removable degenerate solutions for accretion flow 
with multiple critical points, and are shown to be structurally unstable. Depending on suitable 
initial boundary conditions, a homoclinic trajectory can be combined with a standard non 
homoclinic orbit through an energy preserving Rankine-Hugoniot type of stationary shock, 
and multi-critical accretion flow then becomes truly multi-transonic. An effective Lyapunov 
index has been proposed to analytically confirm why certain class of transonic flow can not 
accommodate shock solutions even if it produces multiple critical points. 

Key words: accretion, accretion discs - black hole physics - hydrodynamics - gravitation - 
shock wave - relativity 



1 WHITHER SHOCKED ACCRETION DISC? 

For accretion of matter onto astrophysical black holes, the local radial Mach number Af of the accreting fluid can be defined as the ratio of the 
radial component of the local dynamical flow velocity to that of the propagation of the acoustic perturbation embedded inside the accreting 
matter. The flow will be locally subsonic or supersonic according to Af < 1 or > 1. The flow is transonic if at any moment it crosses the 
Af = 1 hypersurface. This happens when a subsonic to supersonic or supersonic to subsonic transition takes place either continuously or 
discontinuously. The point(s) where such crossing takes place continuously is (are) called sonic point(s), and where such crossing takes place 
discontinuously are called shocks or discontinuities. The particular value of the radial distance r for which M = 1, is referred as the transonic 
point or the sonic point, and will be denoted by hereafter. For r < r^, infalling matter becomes supersonic. Any acoustic perturbation 
created in such a region is destined to be dragged toward the black hole, and can not escape to the domain r > rs. In other words, any 
co-moving observer from r < rs region can not communicate with any observer (co-moving or stationary) located in the sub-domain r > Vs 
by sending any signal which travels with velocity iisignai SjCs, where Cs is defined as the velocity of propagation of the acoustic perturbation 
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(the sound speed) embedded in the moving fluid. Hence the hypersurface through Ts is generated by the acoustic null geodesies, i.e., by 
the phonon trajectories, and is actually an acoustic horizon for stationary configuration, which is produced when accreting fluid makes a 
transition from subsonic (M < 1) to the supersonic {M > 1) state. At a distance far away from the black hole, accreting material almost 
always remains subsonic (except possibly for the supersonic stellar wind fed accretion) since it possesses negligible dynamical flow velocity. 
On the other hand, the flow velocity will approach the velocity of light c while crossing the event horizon, while the maximum possible value 
of sound speed, even for the steepest possible equation of state, would be c/\/3, resulting M > 1 close to the event horizon. In order to 
satis fy such i nner boundary condition imposed by the event horizon, accr e tion onto black hol es exhibit transonic properties in general, see. 



e.g. 



Pringld Jl98lhJchakrabartil Jl990l)jKato. Fukue & Mineshigd 119981) Jprank et al.l j2002h for further details. 



It will be shown in the subsequent sections that the physical transonic accretion solutions can formally be realized as critical solution 
on the phase portrait (spanned by dynamical flow/Mach number and the radial distance) of the black hole accretion. It will also be discussed, 
in great details, that usually a critical point is not identical with a sonic point. For sub-Keplerian angular momentum distribution of matter, 
frequently it happens that the critical features are exhibited more than once in the phase portrait of a stationary solutions describing an 
axisymmetric black hole accretion. For such situations, the number of critical points, unlike spherical accretion, may exceed one. A large 
pool of literat ure, dealing with the theory of the black hole ac cretion disc, studied various properties of such flows wit h mul t iple c r itical points 
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For realistic astrophysical systems, such sub-Keplerian weakly rotating flows are exhibit ed in vari ous physical situations, such as de- 
tached binary systems fed by accretion from OB st ellar winds (Illarionpv & Sunyaev (I975); lLiang & N olan (19^8(D), semi-detached low- 
mass no n-magnetic binari e s ( iBisikalo et al. i I998I) ). and super-massive black holes fed by accretion from slowly rotating central stellar 
clusters jlllarionovl ( Ii988I) ; IHoI jl999l) and references t herein). Even f or a standard Keplerian accretion disc, turbulence may produce such 



low angular momentum flow (see, e.g..llgumenshch ev & Abramowicz 1 1999) , and references therein). 

In supersonic astrophysical flows, perturbation of various kinds may produce discontinuities. Such discontinuities are said to take place 
over one or more surfaces when any dynamical and/or thermodynamic quantity changes discontinuously as such surfaces are crossed. The 
corresponding surfaces are called the surfaces of discontinuity. Certain boundary conditions are to be satisfied across such surfaces, and 
according to those conditions, surfaces of discontinuities are classified into various categories, the most important being the shock waves or 
shocks. Such shock waves are quite often generated in various kinds of supersonic astrophysical flows having intrinsic angular momentum, 
resulting the final state of the flow to be subsonic. This is because the repulsive centrifugal potential barrier experienced by such flows is 
sufficiently strong to brake the infalling motion and a stationary solution could be introduced only through a shock. Rotating, transonic 
astrophysical fluid flows are thus believed to be 'prone' to the shock formation phenomena. 

The non-linear equations describing the steady, inviscid axisymmetric flow in the Kerr metric can be tailored to form a first order 
autonomous dynamical system. Physical transonic solution in such flows can be represented mathematically as critical solutions in the 
velocity (or Mach num ber) phase plane of the flow - they are assoc iated with the critical points (alternatively known as the fixed points or the 
equilibrium points, see lJordan & Smith Jl999h and lChicond for further details about the fixed point analysis techniques). To maintain 

the transonicity such critical points will perforce have to be saddle points, which will enable a solution to pass through themselves. 

Before we proceed further, it is important to introduce a few nomenclatures we will be following in this work. When a configuration 
can have multiple critical points accessible to the accretion flow, we will mention it as 'multi-critical' accretion. Critical points for a non 
dissipative system can be of two types, either saddle or centre. A saddle type critical points allows an accretion flow to pass through it whereas 
a centre type critical point does not. Hence the concept of a 'sonic point' can only be associated with a saddle type critical point and not with 
a centre type one. No transonic solution can pass through a centre type critical point. We will demonstrate that a saddle type critical point 
usually has a sonic point associated with it. As a result, a multi-critical flow having three critical points, two saddle and one centre type, 
can only have two sonic points. A 'multi-transonic' flow, according to our definition, will be a particular subclass of multi-critical transonic 
accretion which can pass through two sonic points. This is possible, as we will see in subsequent sections, only if a shock forms. Otherwise, 
even if a flow has more than one saddle type critical points, transonic accretion passes only through a single sonic point. Such flows will in 
generally be termed as a 'multi-ctritical mono- transonic' flow in our work. On the contrary, when the number of the critical point itself is one, 
and that critical point is of saddle type so that a sonic point is also associated with the flow, we will call it a 'mono-critical mono-transonic' 
flow. We will see in the subsequent sections that among the several phases of the transonic accretion flow, a true multi-transonic flow forms 
only for a shocked accretion, otherwise all accretion flows are mono-transonic: either of the multi-critical mono- transonic category or of 
mono-critical mono- transonic category. All these issues will further be elaborated with greater detail in subsequent sections. 

The true multi-transonicity in the Kerr geometry implies the existence of three critical points and a standing shock, and under the 
fundamental physical requirement of accretion being a process whereby a flow solution should connect infinity to the event horizon of the 
black hole, the three critical points should be such that there would be two saddle points flanking a centre type point between themselves. 
Originating at a large distance, subsonic accretion encounters the outermost saddle type sonic point and becomes supersonic. Subjected to 
the appropriate perturbative environment, such supersonic flow encounters a shock and becomes subsonic again. Such subsonic matter has 
to pass through another saddle type sonic point to meet the aforesaid inner boundary condition. For accretion onto black hole, presence of at 
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least two saddle type sonic points are thus a necessary (but not sufficient) condition for the shock formation. Multi-transonicity, thus, plays a 
crucial role in studying the physics of shock formation and related phenomena in connection to the black hole accretion processes. 

One also expects that a shock formation in black-hole accretion discs might be a general phenomenon because shock waves in rotating 
astrophysical flows potentially provide an important and efficient mechanism for conversion of a significant amount of the gravitational 
energy into radiation by randomizing the directed infall motion of the accreting fluid. Hence, the shocks play an important role in governing 
the overall dynamical and radiative processes taking place in astrophysical fluids and plasma accreting onto black holes. 

The hot and dense post shock flow is considered to be a po werful diagnostic tool in unders tanding various astrophysical phenomena 
like the spectral prope rties of the galactic black hole cand idates dChakrabarti & Titarchuklll995h and that of the supermassive black hole 
at our Ga lactic centre I Moscibrod zka. Das & Czem v 2006), the formation and dynamics of accretion powered galactic and extra-galactic 
outflows JPas & ChakrabartL1999i: 



Das & Chattopadhyaj 2007t) . shock induced nucleosynthesis in the black 
hole accretion disc and the metalicity of the intergalactic matte r (i Mukhopadhy ay 1999) and references therein), and the orig in of the quasi- 
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periodic oscillation in galactic sources (i Sponholz & Molteni 19941 : iDas. Rao & Vadawale 2003 ; Okuda et al. 2004 , 2007 ) and references 
therein). 

The study of steady, standing, stationary shock waves produced in black hole accretion and related phenomena thus acquired an impor- 
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same set of initial boundary conditions describing the multi-transonic accretion, there exist flow profiles which may or may not accommo- 
date shock solutions, the very existence of accretion flow with shock is still a matter of debate to the astrophysical community involved in 
the study of the theory of black hole accretion disc. 

This present treatment confirms, using a much broader and more comprehensive procedure in its scope and objectives than those reported 
previously in the literature, that the shock is an integral part of the black hole axisymmetric accretion in its most generic form, i.e., for the 
general relativistic accretion disc in the Kerr metric. In the present paper we study the sequence of the stationary axisymmetric solutions 
characterized by the angular momentum paying attention to the changes of the flow topology in the extended velocity phase space diagram. 
By the phrase 'extended velocity phase diagram', here we mean the two dimensional graphical representation where the radial Mach number 
has been plotted along the ordinate and the radial distance has been plotted along the abscissa. The radial Mach number is the ratio of 
the radial velocity measured on the equatorial plane of the accretion disc, and the height averaged polytropic sound speed defined on the 
equatorial plane of the disc. The radial distance has been measured from the black hole along the equatorial plane, and is scaled in the units 
of GMbh /c? , where G is the Universal gravitational constant and Mbh is the mass of the black hole considered, G and c have been scaled 
to be unity in the system of geometrical unit. We discuss in detail the critical solution which divides the family of compact transonic flows 
(solutions with sonic point at the innermost critical point) from extended transonic flows (solutions with sonic point at the outermost critical 
point). We argue that the existence of this critical solution is likely to lead to a hysteresis effect when the system chooses the solution with, 
or without a shock, depending on the flow history. In addition, we have proposed a Lyapunov like treatment, to show that the likelihood of 
the event that a shock will form or not, clearly depends on the value of such effective Lyapunov index. 

The plan of the paper is as follows: 

In the next section, we will describe the overall flow configuration, and the corresponding governing equations will be derived, solution 
of which, along with a detailed discussion of the transonic properties arising out of such solutions, will be provided in the section[3]Section|4] 
and|5]will elaborate the multi-transonic behaviour of such flow. In section[6]we will describe certain 'hysteresis effects' corresponding to the 
secular evolution of the flow. Section|7]will illustrate the details of an Lyapunov like treatment for the state transition in transonic accretion 
solution to establish what prompts a multi-transonic accretion to (or not to) contain a standing shock solution. Finally we conclude in the 
section[8] 



2 DRESSING UP THE DISC 

2.1 Spacetime geometry and the stress energy tensor 

To provide a generic description of axisymmetric fluid flow in strong gravity, one needs to solve the equations of motion for the fluid and 
the Einstein equations. As of our present (limited) understanding about the detail physics of the black hole accretion flow is concerned, this 
seems to be an enormously difficult, if not impossible, task to accomplish. The problem may be made relatively tractable by assuming the 
accretion to be non-self gravitating so that the fluid dynamics may be dealt in a metric without back-reactions. The most general form of 
the energy momentum tensor for such non self gravitating compressible hydromagnetic astro physical fluid (with a froze n in magnetic field) 
vulnerable to the shear, bulk viscosity and generalized energy exchange, may be expressed as iNovikov & Thomelll973l) : 



T^-- ^T^J: +T^^- (1) 

where TJJ,'^ and Tg"^ are the fluid (matter) part and the Maxwellian (electromagnetic) part of the energy momentum tensor T'j^ and Tg"^ may 
be expressed as: 
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B'/i"" - 2B''B'') 



(2) 



In the above expression, pv'^v'^ is the total mass energy density excluding the frozen-in magnetic field mass energy density as measured in 
the local rest frame of the baryons (local orthonormal frame, hereafter LRF, in which there is no net baryon flux in any direction), ph^'^ is the 
anisotropic pressure for incompressible gas (had it been the case that 6 would be zero). ? and ri are the co-efficient of bulk viscosity and of 
dynamic viscosity, respectively. Hence ~<;9h'^" and —2i]a^'^ are the bulk viscosity and the viscous shear stress, respectively, q^u" + v'^q" is 
the energy and momentum flux, respectively, in LRF of the baryons. In the expression for Tg", B^/Stt in the first term represents the energy 
density, in the second term represents the magnetic pressure orthogonal to the magnetic field lines, and in third term magnetic tension along 
the field lines (all terms expressed in LRF), respectively. 

Here, the electromagnetic field may be described by the field tensor JT'"' and it's dual jT*'"' (obtained from JT'"' using Levi-Civita 
totally antisymmetric tensor e'"'"'^) satisfying the Maxwell equations through the vanishing of the four-divergence of F*'^" . A complete 
description of flow behaviour could thus be obtained by taking the co-variant derivative of T^" and pv^ to obtain the energy momentum 
conservation equations and the conservation of baryonic mass. 

However, at this stage, the complete solution remains analytically untenable unless we are compelled to adopt a number of simplified 
approximations. Our work concentrates on the inviscid accretion of hydrodynamic fluid. Hence may be described by the standard form 
of the energy momentum (stress-energy) tensor of a perfect fluid: 



T" 



(e +p) UpUy +pgtiv, or, T = (e + p) v®v + pg 



(3) 



Our calculation will thus be focused on the stationary axisymmetric solution of the energy momentum and baryon number conservation 
equations 



= 0; 



[pv" 



= 0, 



Specifying the metric to be stationary and axially symmetric, the two generators = (d/dt)'^ and (j)'^ = (d/d(j))' 
axial isometry, respectively, are Killing vectors. 



(4) 

of the temporal and 



To describe the flow, it is advisable to use the Boyer-Lindquist co-ordinate iBover & Lindquis^ i 1967 ), and an azimuthally Lorentz 
boosted orthonormal tetrad basis co-rotating with the accreting fluid. At this stage we are not interested in non-axisymmetric disc structure, 
hence we neglect any gravo-magneto- viscous non-alignment between the flow angular momentum and the black hole spin angular momen- 
tum. We consider inviscid accretion and denote the constant specific flow angular momentum to be A. Hence, while describing the accretion 
disc dynamics, the viscous transport of the angular momentum has not explicitly been taken into account. Viscosity, however, is quite a 
subtle issue in studying the disc accretion. Even thirty six years after the discovery of standard accretion disc theory dShakura & Sunvaev 
I973I : iNovikov & Thorne|[l973l) . exact modeling of viscous transonic black-hole accretion, including proper heating and cooling mecha- 
nisms, is still quite an arduous task, even for a Newtonian flow, let alone for general relativistic accretion. Nevertheless, extremely large 
radial velocity close to the black hole implies t™/ <C Tvisc, where Tinf and Tvisc are the infall and the viscous time scales, respectively. 
Large radial velocities even at la r ger distances are due to the fact that t he angular momentum co ntent of the accreting fluid is relatively low 
i Beloborodov & Illarionovlll991 : Igumenshchev & Beloborodov 1997 : Proga & Begelmar]|2003 1. Hence, the assumption of inviscid flow is 
not unjustified from an astrophysical point of view. However, one of the most significant effects of the introduction of viscosity would be the 
reduction of the angular momentum. In our work, it has been observed that the location of the sonic points anti-correlates with A, i.e. weakly 
rotating flow makes the dynamical velocity gradient steeper, which indicates that for viscous flow the sonic horizons will be pushed further 
out and the flow would become supersonic at a larger distance for the same set of other initial boundary conditions. Also we found how the 
shock location (and strength, and other related variables) get modified because of the variation of the flow angular momentum. Hence in our 
work, we practically have demonstrated the effect of viscosity as well. 

We consider the flow to be 'advective', i.e. to possess considerable radial three- velocity. The above-mentioned advective velocity, which 
we hereafter denote by u and consider it to be confined on the equatorial plane, is essentially the three-velocity component perpendicular 
to the set of hypersurfaces {Eu} defined by v'^ = const, where v is the magnitude of the 3-velocity. Each Ttv is timelike since its normal 
r/fi oc d^v^ is spacelike and may be normalized as T]''ri^ — 1. 

We then define the specific angular momentum A and the angular velocity as 



A = - 



Vt 



— ^ 



gt4, + Agtt 

g4>4> + M 



40 



The metric on the equatorial plane is given by JNovikov & Thomd ( 1 19730 ) 



ds^ — g^^dx'^dx'^ 



1, z 



-dr 

A 



+ 



- ijjdtf + ^ 
' A 



dr^ + dz^, 



(5) 



(6) 



A 

where A = — 2r + a^, A = r"^ + r^a^ + 2ra^, and lo = 2ar/A, a being the Kerr parameter related to the black-hole spin. The 
normalization condition v'^v^ = —1, together with the expressions for A and Q, in l|5j, provides the relationship between the advective 
velocity u and the temporal component of the four velocity 

-,1/2 

= — im-T^ 7-; , . ^o o o TTTT ■ (7) 



(1 



-4AarA + AV2 (4a2 
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In order to solve we need to specify a realistic equation of state. In this work, we concentrate on polytropic accretion, and some relevant 
details of the flow thermodynamics is presented in the next section. 

It is to be noted here that the use of the polytropic equation of state in describing the thermodynamic properties of the accreting matter 
is quite common in the theory of relativistic black hole accretion. However, one also understands that the polytropi c accretion is not the o nly 
choice to describe such accre tion. Equations o f state other than the adiabatic one, such as the isothermal equation i Yang & Kafatosl[l995 ) or 
the two-temperature plasma l lManmotolEoOoL . have also been used to for this purpose. 



2.2 The thermodynamics of the flow 

The polytropic equation is taken to be of the following form 

P = Kp', (8) 

where the polytropic index 7 (equal to the ratio of the two specific heats Cp and Cv) of the accreting material is assumed to be constant 
throughout the fluid. A more realistic model of the flow would perh aps require a variable polytropic ind ex having a functional dependence 
on the radial distance, i.e. of the functional form 7=7(r); see, e.g., 1 Ryu. Chattopadhvav & Choi 2006i) and references therein for further 
discussions. However, we have performed the calculations for a sufficiently large range of 7 and we believe that all astrophysically relevant 
polytropic indices are covered in our analysis. 

The constant if in l[8j may be related to the specific entropy of the fluid, provided there is no entropy generation in the flow. If in addition 
to (Ull the Clapeyron equation for an ideal gas holds 

P=—PT, (9) 

where T is the locally measured temperature, fi the mean molecular weight, mn^rnp the mass of the hydrogen atom, then the specific 
entropy, i.e. the entropy per particle, is given by ULandau & Lifshitz„1987i) : 

1 7 

a — log K H h constant, (10) 

7 — 1 7 — 1 

where the constant depends on the chemical composition of the accreting material. Equation l llOt confirms that K in l[8j is a measure of the 
specific entropy of the accreting matter 

The specific enthalpy of the accreting matter can now be defined as 

/.= ^, (11) 
P 

where the energy density e includes the rest-mass density and the internal energy and may be written as 

e = P+^-r. (12) 

7-1 

The adiabatic speed of sound is defined by 
2 dp 



c, 

at 



constant entropy 

From M2\ we obtain 

dp fl-l-cl" 



(13) 



ae V 7 ^ 1 
Combination of l ll3t and l[8} gives 

cl = , (15) 

ae 

Using the above relations, one obtains the expression for the specific enthalpy 

h= ^ 7 \ ■ (16) 

The rest-mass density p, the pressure p, the temperature T of the flow and the energy density e may be expressed in terms of the speed of 
sound Cs as 
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T = 



7 



7-1 
1 



7 - 1 - ci 



1 + i 

7 



7-1- 



(19) 
(20) 



We now need to define a specific geometrical structure of the disc, wliich has been done in the next section. 
2.3 The disc geometry 

We assume that the disc has a radius-dependent local thickness H(r), and its c entral plane coinc i des with the equatorial pl a ne of the black 
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Abramowicz et al 



h988l):lchen&Taanil 



1993'):'Kafa tos & Yang|fl994l)jArtemova et alJ ( fl996l)jNaravan. Kato & Honmj jl997l)Jwiital jl998l) ; lHawlev & KroliklhoOlhJArmitage. Reynolds & Chiai 



1 2001) and many others) to use the vertically averaged model in describing the black-hole accretion discs where the equations of motion ap- 
ply to the equatorial plane of the black hole. We follow the same procedure here. The thermodynamic flow variables are averaged over the 
disc height, i.e., a thermodynamic quantity y used in our model is vertically averaged over the disc height as 

__ j"{r)o{ydh) 

We follow I Abramowicz . Lanza & P erviva? to derive an expression for the disc height H{r) in our flow geometry since the relevant 

equations in Abramowicz. Lanza & Perviva l (1997) ars non-singular on the horizon and can accommodate both the axial and a quasi-spherical 
flow geometry. In the Newtonian framework, the disc height in vertical equilibrium is obtained from the z component of the non-relativistic 
Euler equation where all the terms involving velocities and the higher powers of (z/r) are neglected. In the case of a general relativistic disc, 
the vertical pressure gradient in the comoving frame is compensated by the tidal gravitational field. We then obtain the disc height 



H{r) 



7 + 1 



(7-l)c? 



{7-(l + ci)}{A2t;?_a2(^,-l)} 



(21) 



which, by making use of Q, may be be expressed in terms of the advective velocity ^ 



3 GOVERNING EQUATIONS AND THE SOLUTION PROCEDURE 
3.1 The first integrals of motion 



Bhattacharjee & Ray (2007); Goswami et al. (2007); Bhattacharjee et ; 
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From analytical perspective, problems in bl ac k hole accretion fall under the ge neral class of nonl i near dynamicslRav & Bhattac harjee (2002) 
Afshordi & Paczvhski'('2003');'Rav'('2003a'^;'Rav & Bhattacharjee' ('2005a||b| ); | Chaudhurv et alj j2006l) ; lRav & Bhattacharjee! ( I2006i.i2007an 



since accretion describes the dynamics of a compressible 



astrophysical fluid, governed by a set of nonlinear differential equations. Physical transonic solutions can be mathematically realized as 
critical solutions in the phase portrait of the flow. To obtain such critical solutions, it is first convenient to construct a set of first integrals of 
motion. 

The temporal component of the energy momentum tensor conservation equation leads to the constancy of the total specific energy of 
the accretion flow along each streamline of the flow. This specific energy is denoted by £ (which actually is t he rel ativistic analogue of the 
Bernoulli's constant, and can be expressed in terms of the flow enthalpy and vt, as £ = hvt, see lAndersonI Jl989l) for further detail), and 
hence from Q and J16b it follows that: 



£ 



(7- 



7-(l + c2) 



Ar^A 



The rest-mass accretion rate M is obtained by integrating the relativistic continuity equation. One finds 
M = 4-kA^Hp- ^ 



Here, we adopt the sign convention that a positive u corresponds to accretion. The entropy accretion rate H can be expressed as: 



AtvA 2 c. 



VT 



(7- 



7-(l + ci) 



Hir) 



(22) 



(23) 



(24) 



One can solve the conservation equations for £ , Al and H to obtain the complete accretion profile. 

We thus have two primary first integrals of motion along the streamline - the specific energy of the flow £ and the mass accretion rate 
M. Even in the absence of creation or annihilation of matter, the entropy accretion rate H is not a generic first integral of motion. As the 
expression for H contains the quantity K=p/ , which is a measure of the specific entropy of the flow, the entropy accretion rate H remains 
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constant throughout the flow only if the entropy per particle remains locally invariant. This condition may be violated if the accretion is 
accompanied by a shock. Thus H is conserved for shock free polytropic accretion (and wind) and becomes discontinuous (actually, increases) 
at the shock location, if such a shock is formed. However, H is of utmost importance for our purpose, since it acts as the degeneracy remover 
for the multi-critical accretion, see subsequent discussions for further detail. 

Being equipped with the disc geometry and the integrals of motion, we will now discuss the transonicity of the flow in the next section. 



3.2 Velocity gradient and the transonicity 

The gradient of the acoustic velocity can be computed by differentiating ( 124b and can be obtained as: 



dr 
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The dynamical velocity gradient can then be calculated by differentiating ilH with the help of ( 1251 ) as: 
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(26) 



where (see dBarai et al ]|2004h for further detail): 



= \ Vt — a {vt — I) , Ipa = [I 
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+ \gt4,) 

A real physical transonic flow must be smooth everywhere, except possibly at a shock. Hence, if the denominator of the equation ( |26| ( 
vanishes at a point (for certain value of r), the numerator must also vanish at that point to ensure the physical continuity of the flow. One 
thus arrives at a critical point (or the fixed point/equilibrium point - by borrowing the terminology used in the theory of dynamical systems) 
conditions by simultaneously making the numerator and the denominator equal to zero. The critical point conditions can thus be obtained as: 



u^ (7 + 1) 



2il] — u^VtCJ 



1/2 



■I(r=re) 



xAr 



2r (r - 1) + 4A 
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(28) 



For any value of ,A,7,a], substitution of the values of it|(r=rc) ^rid Cajr=r^ in terms of rc in the expression for £ M2\ . provides a 
polynomial in r^ the solution of which determines the location of the critical point(s) rc- 



3.3 Non isomorphism of the critical and the sonic points 



It is obvious from ( |28t that u\^^_^ ^ -^Csr^v^^ and hence the Mach number at the critical point is jio/ equal to unity in general. This phenomena 
can more explicitly be demonstrated for relativistic disc accretion onto a Schwarzschild black hole. For Schwarzschild metric, we calculate 
the Mach number of the flow at the critical point as 



where 



h{rc.\) 



7 + iy h{rc,\) + h{rc,\) 



3ri - 2Vrc + 3V 



f2(rc,\) = 



2rc 



2ri 



AVe + A^ 



(29) 



(30) 



_ - X^rdrc - 2) ' ■' ' ' ' rc(rc-2) ri ~ \^rc{ra- 2) 

Clearly, A/c is generally not equal to unity, and for 7 ^ 1, is always less than one, and it can easily be shown that for a realistic choice of 
initial boundary conditions, the separation between the critical and the sonic points can be as high as several hundred gravitational radii. . 

Hence, the sonic points are not isomorphic to the critical points in general, neither numerically, nor topologically, and we categorically 
distinguish a sonic point from a critical point. In the literature on transonic black-hole accretion discs, the concepts of critical and sonic points 
are often made ident i cal by erroneously definin g an 'effective' sound speed leading to the 'effective' Mach number (for further details, see, 
eg. 



Matsumoto etall d 1 9841) : [chakrabartil d 1 989l) ) for an 'effective' geometry. For realistic flow in general relativity, the exact sonic condition 



should be obtained for at the location where the bulk flow velocity has to match the local unsealed speed of propagation of the acoustic 
perturbation. We thus strongly disagree to accept the synonymous use of the terms critical and sonic points for the following reasons: 

In the existing literature on non general relativistic pseudo-Schwarzschild transonic disc accretion, the Mach number at the critical poin t 
turns out to be a function of 7 only, and hence Mc remains constant if 7 is constant. For example, use of the IPaczvnski & Wiitj jl980l) 
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pseudo-Schwarzschild potential to describe tlie adiabatic accretion phenomena leads to (see iDasI ( 120020 ■ and references therein for further 
details of such calculation) 



The above expression does not depend on the location of the critical point and depends only on the value of the adiabatic index chosen to 
describe the flow. Note that for isothermal accretion 7 = 1, and hence the sonic points and the critical points are identical. 

However, the quantity in Eq. l |29l ) as well as when calculated for ( 1281 ), is clearly a function of and hence, generally, it takes 
different values for different for transonic accretion. The Mach numbers at the critical (saddle type) points are Ihus functional , and not 
functions, of the initial boundary conditions. The difference between the radii of the critical point and the sonic point may be quite significant. 
One defines the radial difference of the critical and the sonic point (where the Mach number is exactly equal to unity) as 

Arc = \rs — rc\. (32) 

The quantity Ar^ may be a complicated function of [£, X, 7, a], the form of which can not be expressed analytically, but it can easily be 
evaluated using numerical integrations, and we have done so. We estimate the Ar" can be as large as 10^ rg, or even more, for flow passing 
through the outer critical/sonic points. 

The radius rs in Eq. B2t is the radius of the sonic point corresponding to the same , A, 7, a] for which the radius of the critical point 
rc is evaluated. Note, however, that since rs is calculated by integrating the flow from rc upto the point where the Mach number exactly 
becomes unity, Ar^ is defined only for saddle-type critical points. A physically acceptable transonic solution can be constructed only through 
a saddle-type critical point, and not through a centre type critical point. Hence the concept of a sonic point corresponding to a centre type 
critical point is meaningless, and multi-transonic accretion possesses two sonic points (and three critical points). 

Lemma 3.1. Physically acceptable global transonic solutions must always produce odd number of critical points, whereas even number of 
sonic points are produced for the multi-transonic accretion. Odd number of sonic point is observed for mono-transonic accretion, however, 
this is a trivial case because the number of sonic point is unity. 

However, for saddle-type critical points, rc and rs should always have one-to-one correspondence, in the sense that every critical point 
that allows a steady solution to pass through it is accompanied by a sonic point, generally located at a different radial distance. 

It is worth emphasizing that the distinction between critical and sonic points is a direct manifestation of the non-trivial functional 
dependence of the disc thickness on the fluid velocity, the sound speed and the radial distance, i.e., on the disc geometry as well as the 
equation of state in general. In the simplest idealized case when the disc thickness is assumed to be constant, or the infalling matter being 
described using the isothermal equation of state, one would expect no dist inction betwee n critical and sonic points. In this case, as has 
been demonstrated for a thin disc accretion onto the Kerr black hole .Abraham et al. ( 20061) . the quantity Arc vanishes identically for any 
astrophysically relevant value of [£, A, 7, a]. 



4 MARCH TOWARD THE MULTI-TRANSONICITY 
4.1 The parameter space diagram 

By substituting the value of the bulk velocity and the acoustic velocity evaluated at the critical points to the expression for the energy first 
integral, one obtains a quasi polynomial (not an exact polynomial in strict mathematical sense, since it contains fractional power terms) in rc 
parametrized by [£, A, 7, a]. Solution of such a quasi polynomial (with additional physical constraint that all solutions must be real, positive 
and will be greater than rg), provides the critical points for general relativistic axisymmetric accretion in the Kerr metric. One finds at most 
three critical points for relativistic disc accretion for some values of [£, A, 7, a]. 

The above mentioned quasi polynomial, and hence the critical points, are completely determined by four parameters [£,X,j,a]. A 
four dimensional parameter hyper-space is thus required to be looked upon to obtain the global set of transonic solutions. For the shake of 
convenience, we analyze a two dimensional projection of such four dimensional hyperspace. Since £, A, 7, a are mutually orthogonal (i.e., 
choice of one among those four parameters does not influence the choice of the other parameter(s)), a total number of *C2 different choices 
are allowed for selecting such projections, each characterized with two different parameters out of the four parameter set , A,7,a], by 
keeping the other two parameters at a pre defined fixed value. In this work, we choose to project the hyperspace on an [£ — A] plane by 
keeping 7 = 4/3 and a = 0.3. However, similar projected [£ — A] submanifold can routinely be analyzed for other values of 7 and a. 

To begin with, we first set the astrophysically relevant bounds on {£,"/} to model the realistic situations encountered in astrophysics. 
Since the specific energy £ includes the rest-mass energy, £ = 1 is the lower bound which corresponds to a flow with zero thermal energy 
at infinity. Hence, the values £ < 1, corresponding to the negative energy accretion states, would be allowed if a mechanism for a radiative 
extraction of the rest-mass energy existed. The possibility of such an extraction would in turn imply the presence of the viscosity or other 
dissipative mechanisms in the fluid. Since we concentrate only on non-dissipative flows, we exclude £ < 1. On the other hand, although 
almost all f > 1 are theoretically allowed, large values of £ represent flows starting from infinity with very high thermal energy. In particular. 
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Figure 1. The two dimensional parameter space division for low angular momentum axisymmetric accretion in the KeiT metiic. The parameter space is 
spanned by the conserved specific energy of the flow £ and the specific angular momentum A, and has been drawn for a fixed value of 7 = 4/3 and a = 0.3. 
The O and the I regions represents the mono-critical mono-transonic accretion solutions passing through the saddle type outer and the inner critical (sonic) 
points, respectively. The region marked by A represents multi-critical accretion, and can further be divided (not shown in the figure) into two different regions, 
A|\is where shock does not form even if the flow has three ciitical points, and Ag, which accommodates true multi-transonic accretion with general relativistic 
Rankine Hugoniot shock. The region marked by Ai represents mono-transonic accretion passing through the inner sonic point and accompanied by an 
additional pair of non-accessible critical points and a homoclinic orbit through one of those, see text for further details. 

£ > 2 accretion represents enormously hot flow configurations at very large distance from the black hole, which are not properly conceivable 
in realistic astrophysical situations. Hence, we set 1 ^ f 2. 

The physical lower bound on the polytropic index is 7 — 1, which corresponds to isothermal accretion where accreting fluid remains 
optically thin. Hence, the values 7 < 1 are not realistic in accretion astrophysics. On the other hand, 7 > 2 is possible only for superdense 
matter with a very large magnetic field and a direction-dependent anisotropic pressure. The presence of a magnetic field would in turn require 
solving the general relativistic magneto-hydrodynamic equations, which is beyond the scope of this paper. Thus, we set 1 7 ^ 2. However, 
astrophy sically preferred values of 7 for realistic black-hole accretion range from 4/3 (ultra-relativistic) to 5/3 (purely non-relativistic flow) 
I Frank e t al. 2002). Hence, we mainly focus on the parameter range 



In figure 1, the regions marked by O and I correspond to the formation of a single critical point, and hence the mono-transonic disc 
accretion is produced for such region. In the region marked by I, the critical points are called 'inner type' critical points since these points 
are formed sufficiently close to the event horizon, in most of the cases even closer than the innermost stable circular orbit (ISCO). In the 
region marked by O, the critical points are called 'outer type' critical points, because these points are located relatively far from the black 
hole. Depending on the value of [£,X, 7, a], an outer critical point may be as far as 10*' rg, or more. 

The outer type critical points for the mono-transonic region are formed, as is obvious from the figure, for sufficiently weakly-rotating 
flows. For sufficiently low angular momentum, accretion flow contains less amount of rotational energy, thus most of the kinetic energy 
in utilized to increase the radial dynamical velocity u at a faster rate, leading to a higher value of du/dr. Under such circumstances, the 
dynamical velocity u becomes large enough to overcome the acoustic velocity Cs at a larger radial distance from the event horizon, leading 
to the generation of supersonic flow at a large value of r, which results the formation of the sonic point (and hence the corresponding critical 
point) far away from the black hole event horizon. On the contrary, the inner type critical points are formed, as is observed from the figure, 
for strongly rotating flow in general. Owing to the fact that such flow would posses a large amount of rotational energy, only a small fraction 
of the total specific energy of the flow will be spent to increase the radial dynamical velocity u. Hence for such flow, u can overcome Cs only 
at a very small distance (very close to the event horizon) where the intensity of the gravitational field becomes enormously large, producing 
a very high value of the linear kinetic energy of the flow (high u), over shedding the contribution to the total specific energy from all other 
sources. However, from the figure it is also observed that the inner type sonic points are formed also for moderately low values of the angular 




(33) 
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momentum as well (especially in the region close to the vertex of the wedge shaped zone marked by Ai). For such regions, the total conserved 
specific energy is quite high. In the asymptotic limit, the expression for the total specific energy is governed by the Newtonian construct, and 
one can have: 

£-(y) +(^) +('fW..o„. (34) 

\ / linear \ > / tncrmal \ / rotational 

where <& is the gravitational potential energy in the asymptotic limit. From l l34b it is obvious that at a considerably large distance from the 
black hole, the contribution to the total energy of the flow comes mainly (rather entirely) from the thermal energy. A high value of £ (flow 
energy in excess to its rest mass energy) corresponds to a 'hot' flow starting from infinity. Hence the acoustic velocity corresponding to the 
'hot' flow obeying such outer boundary condition would be quite large. For such accretion, flow has to travel a large distance subsonically 
and can acquire a supersonic dynamical velocity u only at a very close proximity to the event horizon, where the gravitational pull would be 
enormously strong. 

The [£, A] corresponding to the wedge shaped regions marked by A and Ai produces three critical points, among which the largest and 
the smallest values correspond to the saddle type, the outer r°"' and the irmer r^", critical points respectively. The centre type middle critical 
point, r^*'', which is unphysical in the sense that no steady transonic solution passes through it, lies in between and r°"'. 

For the kind of accretion flow considered in this work, the critical points of the phase trajectories can be identified first, following 
which a linearized study in the neighbourhood of these critical points may be carried out, to develop a complete and rigorous mathematical 
classification scheme to identify whether a critical point is of saddle type or of centre type. Global understanding of the flow topologies will 
then necessitate a full numerical integration of the non-linear equations of the flow, which has successfully been performed in this work. For 
the shake of completeness and to have a better understanding of the multi-transonic behaviour, it will not be unjustified to briefly discuss the 
classificat ion scheme for the var ious kind of critical points which may be produced in relativistic accretion disc around a Kerr black hole. 



FoUowing lGoswami et al.l 120071), the methodology for performing such classification is presented below. 



4.2 Classification scheme for the critical points 

Eq. ( I26t could be reformulated as 
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with the primes representing the total derivatives with respect to r, and r is an arbitrary mathematical parameter. Here, 

. Ar^A 
fir) 



(35) 



(36) 
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The critical conditions are obtained with the simultaneous vanishing of the right hand side, and the coefficient of d{u^)/dr in the left 
hand side in i35\ . This will provide 
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7 + 1 1*2 

as the two critical point conditions. Some simple algebraic manipulations shows that 

nl = (39) 
f9i 

following which c2|r=r^ can be rendered as a function of rc only, and further, by use of (122b.. rc, c^c and '^c can all be fixed in terms of the 
constants of motion like £, 7, A and a. Having fixed the critical points it should now be necessary to study their nature in their phase portrait 
of versus r. To that end one applies a perturbation about the fixed point values, going as. 



= u^\i=ia + Su^, Cs = c^U^Ta + Sc^s, r = rc + Sr 
in the parametrized set of autonomous first-order differential equations, 

df -7+1191 92 \ dr j\ f 
and 



(40) 



(41) 
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with r being an arbitrary parameter. In the two equations above Sc?^ can be closed in terms of 5v? and 5r with the help of l |25t . Having done 
so, one could then make use of solutions of the form, 5r ~ exp(nr) and 5v? ~ exp(nr), from which, Q, would give the eigenvalues — 
growth rates of 5v? and 5r in f space — of the stability matrix implied by J4U42t . Detailed calculations will show the eigenvalues to be 

= |^'4x?+66L_, (43) 

where = -::p^ and Xi,£,i and ^2 can be expressed as polynomials of Tc. Cl^ can be evaluated for any [E, A, 7, ajonce the value of the 
corresponding critical point rc is known. The structure of ( 143 1 immediately shows that the only admissible critical points in the conserved 
Kerr system will be either saddle points or centre type points. For a saddle point, Vl^ > 0, while for a centre-type point, Q,^ < 0. 

For multi-critical flow characterized by a specific set of [£, A, 7, a], one can obtain the value of Q,^ to be positive for the inner and the 
outer critical points, showing that those critical points are of saddle type in nature. CP comes out to be negative for the middle critical point, 
confirming that the middle critical point is of centre type and hence no transonic solution passes through it. One can also confirm that all 
mono-transonic flow (flow with a single critical point characterized by [£ , X] taken either from I or from O region) corresponds to saddle 
type critical point. 



4.3 Distinction between the two different category of flows, both having three distinct critical points 

There are distinct topological differences between the multi-transonic flow characterized by [£, A] taken from the region marked by A, and 
the region marked by Ai. For region marked by A, the entropy accretion rate H for flows passing through the inner critical point is greater 
than that of the outer critical point 

H (rr) > H (r-r*) (44) 
while for the region marked by Ai, the following relation holds 

H (rr) < H {rT) (45) 

The above two relations show that [£, A] region marked by A represents the possibility of having multi-critical multi-transonic accretion, 
while [£, X]£ , A]^^ corresponds to the mono-transonic accretion, with an additional pair of critical points. For the multi-critical accretion, 
as we will see in more detail in the subsequent sections, and as is illustrated in the panel diagrams (ii) - (iv) of figure 2, one can have a 
possibility that the global accretion solution may pass through the outer critical/sonic point, and a local accretion solution passing through 
the inner critical/sonic point, which folds back onto itself (at a point of inflexion formed at a certain radial distance r=r {£,\, 7, a)), forming 
a homoclinic orbi^ passing through the saddle type inner critical point and embracing the centre type middle critical point flanked by the 
saddle type outer and the inner critical points. Such a homoclinic orbit can be connected with the standard non-homoclinic path through a 
standing shock, and hence, for a suitable initial boundary condition determined by [£,\, 7, a], accretion flow can pass through both the inner 
and the outer saddle type critical/sonic points. Hence we designate such flows (characterized by the set of parameters taken from the region 
A of the parameter space, as shown in figure 1), which can have more than one sonic point available to it to pass through (in reality whether 
it will pass through both the sonic points depends on the possibility of the shock formation, see subsequent sections for further details) to be 
'multi-transonic' accretion. 

The region A, thus, as we will further discuss in the subsequent sections in even more details, can further be divided into two regions. 
For the region where three critical points exist but no shock solution is allowed, the accretion is of the category multi-critical mono-transonic 
in strict sense, and such a region will hereafter be marked as Ams - which reads 'multi-critical mono-transonic accretion flow with no shock'. 
For the rest of the A region, a steady standing Rankine Hugoniot kind of shock forms, and accretion passes through both of the outer and 
the inner sonic points. Accretion is thus truly multi-transonic, and such a region in parameter space will be marked as As, which reads 
'multi-critical multi-transonic accretion with shock'. 

On the other hand, flows described by the parameters taken from the region Ai of the parameter space behaves in a completely different 
way. Here, the global accretion solution passes through the inner critical/sonic point and hence it is mono-transonic. The homoclinic orbit 
is formed through the outer critical point, and the accretion sector of such homoclinic orbit can not be connected with the corresponding 
standard transonic accretion solution passing through the inner sonic point, by a shock (or through by whatever physical means), since 
the corresponding segment (lying in the length scale spanning the region between the outer critical point and the point of inflexion of the 
homoclinic loop) of the global transonic solution remains subsonic (shock can not form in a subsonic flow). We thus call such flow topologies 
as the 'mono-transonic accretion with an additional pair of critical points'. 

It is to be noted that the above mentioned solutions are commonly known as the 'multi-transonic wind' in the literature, owing to the 
fact that for a limited subset of values of [£, A, 7, a] taken from the region Ai, the wind segment of the homoclinic orbit can be connected to 



^ In the termi nology of the autonomo u s dynamical sys tems, a homoclinic orbit on a phase portrait is a solution that connects a saddle type critical point to 
itself. See, e.g.. |jordan & SmitlJ JT999h ; IChiconeH2006l) and the references therein for further detail. 
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the non homoclinic wind like solutions (passing thirough the inner sonic point) through a standing shock. However, we would not like to use 
such terminology here. Such terminology, as we feel, is misleading in the sense that the wind solution have no role in studying the accretion 
phenomena (except, perhaps, for studying the accretion powered outflow), and hence the term 'multi-transonic wind' sounds as if there is no 
accretion taking place for the initial boundary conditions chose from the region Ai, and such nomenclature is of compromised clarity only. 

The actual situation, in contrary, is that there are several phases of transonic accretion flow, as we see in the subsequent section, and we 
concentrate only on the accretion part when we study accretion onto the black holes. The only true multi-transonic accretion are obtained if 
and only if shocks are formed. Otherwise, the accretion is always mono-transonic, whether the flow possesses more than one critical points or 
not. Either such mono-transonic solution has only one saddle type critical point, and one corresponding sonic point (outer or inner), or it can 
have two sonic points available, and a homoclinic orbit passing through the inner critical point, but due to the absence of the shock formation 
phenomena, the transonic flow is realizable only through the outer sonic point, and the flow is essentially topologically isomorphic to the 
mono-critical mono- transonic accretion from the region marked by O in the parameter space as shown in the figure 1, or, in its last variant, 
the flow can have three critical points, but the transonic accretion passes through only through the inner critical/sonic point, accompanied 
by a homoclinic loop through the outer critical point, and the corresponding accretion flow configuration is isomorphic to the mono-critical 
mono-transonic accretion passing through the inner critical/sonic point, as observed for flows chosen from the region marked by I in the 
parameter space as shown in the figure 1. All such issues will be further clarified in greater detail in the next section. 

The boundary between the region marked by A and Ai is a very special one, since for [£, A] defined by the boundary, the entropy 
accretion rate corresponding to the solution passing through the outer critical point is exactly same as that of through the inner critical 
point (H (r^") = H (r°"')). [£, A] defined by this boundary (rather \£, A, 7] for a fixed a, defined by this surface in general) provides the 
non-removable degenerate bistable/unstable solutions, we will discuss this issue in more details later. 

There are other regions for [S , A] space for which either no critical points are formed, or two critical points are formed. These regions 
are not shown in the figure, since none of these regions is of our interest. If no critical point is found, it is obvious that transonic accretion 
does not take place for those set of [£ , A] . For two critical point region, one of the critical points are always of centre type, since according 
to the standard dynamical systems theory two successive critical points can not be of same type (both saddle, or both centre). Hence the 
solution which passes through the saddle type critical point would encompass the centre type critical point by forming a loop like homoclinic 
orbit, and hence such solution would not be physically acceptable accretion solutions since such solutions do not connect infinity to the event 
horizon. In subsequent sections, we will thoroughly describe the methodology for obtaining such flows parametrized by [f, A,7,a] taken 
from the regions A and Ai. 



5 EFFECT OF CHANGE OF THE CONTROL PARAMETER ON THE PHASE PORTRAIT 

5.1 Mono-transonic flow: Characteristics of the integral curves and the methodology to obtain the flow topology 

As has already been mentioned, we choose an [£ — A] projection of the entire four dimensional parameter space. The transonic flow properties 
and the corresponding flow topologies depends on four parameters of the [£, A, 7, a]set. We would like to illustrate such dependence by 
varying one parameter while keeping the other three parameters fixed. We will smoothly vary the angular momentum A of the flow and will 
observe its consequence on the phase topology. The corresponding topologies have been shown in six different panels (i) - (vi) in figure 2. 
All the panel figures are drawn by plotting the Mach number along the Y axis and the distance (measured from the event horizon and scaled in 
the unit of the gravitational radius) along the X axis in logarithmic scale. For all the topologies in the figure, the value of £ has been chosen to 
be 1.0000033, which is in accordance with t he value of the Bernoulli's constant for accretion flow onto the supermassive black hole located 
at our Galactic centre Moscibro dzka. Das & Czemv hood) . The value of 7 has routinely been taken to be 4/3, and a = 0.3 has been taken 



as a representative value of the Kerr parameter. Same procedure can be performed using any other value of a as well. Stable mono-transonic 
global solutions through the outer critical point is available for < A^Ai, where Ai has the value in between 2.839 and 2.84. For such an 
interval of A, the solution topology has been shown in the panel (i) of the figure 2, hereafter described as panel (ii). 

We now describe (in somewhat great detail) the methodology for obtaining a transonic solution topology on a phase portrait. This 
requires numerical integration of the equations describing the velocity gradient, as well as the gradient of the acoustic velocity, of the 
accretion flow. One thus needs to define the value of both the velocity gradients at the critical point, pick up the required quantities (the value 
of the dynamical and the acoustic velocities, and their respective gradients at the critical point, and of the critical point itself), and to use 
those values as the starting values for the integration, like what is to be done in case of the standard initial value problem corresponding to 
the numerical solution of the differential equation. 

To obtain the dynamical velocity gradient at the critical point, one applies I'Hospital's rul e on j26b. After som e algebraic manipulations. 



the following quadratic equation is formed, which can be solved to obtain {du/dr)^^^_^ ^ (see lBarai et al.l <2004l) for further details) 
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Figure 2. Solution topologies for various cases as obtained by gradually increasing the flow angular momentum, by keeeping the other initial boundary 
conditions at a fixed value of [S = 1.0000033, 7 = 4/3, a = 0.3]. The usual accretion branch is denoted by solid line (red coloured in the online version) 
marked by 'Acc', and the other solid line intersecting the accretion branch (at the outer critical points for panel (ii) - (iv) and at the inner critical point at panel 
(vi)) represents the corresponding wind solution. The homoclinic orbit passing through the inner (for panel (ii) - (iv) ) or the outer (for panel (vi)) critical point 
is shown by the dotted line (blue coloured in the online version). The shock transition (at panel (iv) for true multi-transonic accretion, and at panel (vi) for 
shock formation in wind solution) is shown by vertical long dashed Une (green coloured in the onUne version). 
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Note that all the above quantities are evaluated at the critical point. 
Hence we compute the critical advective velocity gradient as 
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where the '+' sign corresponds to the accretion solution and the '-' sign corresponds to the wind solution, see the following discussion for 
further details. Similarly, the space gradient of the acoustic velocity dcs/dr and its value at the critical point has also been calculated. 
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For the solution topologies shown in the figure 2, using the specified set of [£,A,7,a], we first solve the equation for the energy 
first integral defined at the critical point, to obtain the corresponding critical point Vc = 99851.8623, which is the point of intersection in 
between the accretion branch marked by 'Acc' in the figure, and the corresponding wind branch, the other curve shown in the figure. We then 
calculate the critical value of the advective velocity gradient at rc from l l48b . By integrating the equation for the dynamical and the acoustic 
flow velocity gradient (Eq. l l25b and ( 126b . respectively), from the critical point, using the fourth-order Runge-Kutta method, we then calculate 
the local advective velocity, the polytropic sound speed, the Mach number, the fluid density, the disc height, the bulk temperature of the flow, 
and any other relevant dynamical and thermodynamic quantity characterizing the flow. In this way we obtain the accretion branch 'Acc' by 
employing the above mentioned procedure. 

The solution, as is obvious from the figure, is two-fold degenerate owing to the ±it degeneracy which reflects the physical accretion/wind 
degeneracy. We have, however, removed the degeneracy by orienting the curves, and thus each line represents either the wind or accretion. 
We have arbitrarily assigned the + sign solution in Eq. ( 148b to the accretion and the — sign solution in Eq. J48b to the 'wind' branch. This 
wind branch is just a mathematical counterpart of the accretion solution (velocity reversal symmetry of accretion), owing to the presence of 
the quadratic term of the dynamical velocity in the equation governing the energy momentum conservation. 

The term 'wind solution' has a historical origin. The solar wind sol ution first intr oduced by Parker i ParkeJ 1965 ) has the same topology 
profile as that of the wind solution obtained in classical Bondi accretion i BondJl952h . Hence the name 'wind solution' has been adopted in a 
more general sense. The wind solution thus represents a hypothetical process, in which, instead of starting from infinity and heading towards 
the black hole, the flow generated near the black-hole event horizon would fly away from the black hole towards infinity. The topology of 
such a process is represented by the wind solution in the panel, which, as already been mentioned, is the curve intersecting the accretion 
branch 'Acc' at the outer critical point. 

The above procedure for obtaining the flow topology is also applied to draw the mono-transonic accretion/wind branch through the inner 
sonic point. In fact, the same procedure may be used to draw real physical transonic accretion/wind solutions passing through any acceptable 
saddle-type critical point. Note, however, that the sector of the accretion branch starting from the intersection of the accretion and the wind 
solution in the panel diagram and ending at the event horizon is not the complete subsonic branch, because the point of intersection is a 
critical point and not a sonic point. Using the procedure described above we have to integrate the flow from the critical point to the sonic 
point where the Mach number becomes unity. 



5.2 Appearance of the homoclinic orbit and the emergence of the multi-transonicity 

We now further increase A so as to reach a value A2 - greater than but infinitesimally close to Ai, i.e.: 

Ltt^o IA2 — All = e 

A2 thus characterizes the flow topology generated at the right side of the boundary between the O and the A regions shown in figure I, and 
infinitesimally close to it. Such solution has been shown in the panel (ii). The line separating the O and the A region is actually the line 
of bifurcation of the critical point, because the number of critical point increases after crossing such line. The line of separation between 
the O - A region, in general, is a three dimensional hypersurface for the generic multi dimensional parameter space spanned by [£, A, 7, a]. 
Additional critical points appear in pairs once such a line is crossed from the left side (by gradually increasing the value of the flow angular 
momentum A) of the parameter space. Hence, while increasing A smoothly, initially there is a region of a single saddle point, then followed 
by the pair creation of a saddle and a centre type point after crossing the O - A boundary. That is what is manifested in panel (ii). In addition 
to a transonic accretion solution passing through the outer critical point, which is denoted by the solid lines (red coloured in the online 
version) a tiny loop of homoclinic orbit (shown using dotted line, blue coloured in the online version) appears representing the flow through 
(and around) a nascent saddle-centre pair of critical points. Physically, this is what it has to be, because with a centre type point but without 
another saddle point, the flow solutions will all curl about the centre type point, and there will be no means of connecting the event horizon 
with infinity through an isolated solution passing through the saddle type inner critical point, and embracing the centre type middle critical 
point (here, the inner and the middle critical points are formed at 7.2 and 7.4 r^, respectively). The homoclinic path has its existence in 
isolation, and unlike the transonic solution passing through the outer critical point, can not be regarded as a global transonic solution on its 
own. A physically acceptable flow through the inner critical point has to be connected with the solution passing through the outer critical 
point. As we will see in the subsequent paragraphs, a standing Rankine-Hugoniot shock perfectly accomplish the task. 

Accretion flow through the inner critical point has higher value of the entropy accretion rate compared to that of passing through the 
outer critical point. The ratio of Hi„ and Hout attains its maximum value close to the O - A boundary, and then monotonically decreases 
as one smoothly transits to the higher end of A (toward the boundary separating the A and the Ai regions, respectively). 

As we further process toward higher value of A, the area enclosed within the homoclinic orbit starts increasing and finally it fills up the 
entire region bounded by the transonic accretion and wind passing through the outer critical point. The situations have been demonstrated 
sequentially in the successive panel figures (iii), (iv) and (v). 



Hysteresis in shocked black hole accretion 15 



5.3 Multi-critical accretion with and without shock 

In panel (iii), R^. is 212.52. Shock, however, does not form even if three critical points are present. Such multi-critical flow (which is actually 
a mono-transonic accretion because in absence of the shock, accretion flow can pass only through the outer sonic point even if the inner sonic 
point formally exists) without a standing shock is observed for a range of value of A^Aa where A3~3.239, which defines the boundary 
of the Ans region. 

As A is further increased so that A > A3, Rankine-Hugoniot conditions (relativistic Rankin-Hugoniot conditions are explicitly derived 
in the next section) gets satisfied, and standing shock forms. Panel (iv) shows the topology of the multi-transonic accretion with shock 
solutions, where the shock transition has been marked by the vertical long-dashed (green coloured in the online version) line segment. The 
shock location, shock strength, in which is defined as the pre (-) to the post (-I-) shock ratio (defined not in the temporal sequence but in 
the spatial sequence, of the Mach number) and is denoted by M-/M+, and the entropy enhancement ratio at the shock Rj- are found to be 
18.3336735 and 74.9621923, respectively. 

It is to be noted that the panel (iv) is not the solution corresponding to the minimum value of A for which the shock forms, i.e., does not 
correspond to the value of A separating the Ans and the As region. The strongest shock (forms for the maximally allowed value of R= in 
As) are formed for A4 where 

LtE^o|A4 — A3 1 = e 

For example, for A = 3.24, the shock forms at 27.5rg, and the corresponding shock strength and the entropy enhancement ratio at the shock 
are found to be 41.7 and 197.62, respectively. 

The relativistic shock equations and its derivation, along with a detail description of the methodology of finding the multi-transonic flow 
with and without shock is presented in the subsequent paragraphs. However, from panel (iii) and (iv), it is to be understood that shock does 
not necessarily form for all multi-critical flow. Having two saddle type critical points is thus a necessary but not a sufficient condition. The 
reason behind this disparity will be analytically justified in the section ^ and section Q 

5.4 Methodology for obtaining multi-transonic topology 

Using a specific set of [£, A, 7, a], one first solves the equation for E at the critical point to find out the corresponding three critical points, 
saddle type inner, centre type middle, and the saddle type outer. The space gradient of the flow velocity as well as the acoustic velocity at 
the saddle type critical point is then obtained. Such u\^^_^ ^ , Cs\^^_^ ^,dcs/dr\^^_^ ^ and du/dr\^^_^ ^ serve as the initial value condition 
for performing the numerical integration of the advective velocity gradient using the fourth-order Runge-Kutta method. Such integration 
provides the outer sonic point (located closer to the black hole compared to the outer critical point, since the Mach number at the outer 
critical point is less than unity), the local advective velocity, the polytropic sound speed, the Mach number, the fluid density, the disc height, 
the bulk temperature of the flow, and any other relevant dynamical and thermodynamic quantity characterizing the flow. The corresponding 
wind solution through the outer sonic point is obtained in the same way, by taking the other value (note that being a quadratic algebraic 
equation, solution of J48t provides two initial values, one for the accretion and the other for the wind) of du/dr\^^_^ ^ . 

The respective accretion and the wind solutions passing through the inner critical point are obtained following exactly the same proce- 
dure as has been used to draw the accretion and wind topologies passing through the outer critical point. Note, however, that the accretion 
solution through the inner critical point folds back onto the wind solution which is a homoclinic orbit through the inner critical point encom- 
passing the centre type middle critical point. A physically acceptable transonic solution must be globally consistent, i.e. it must connect the 
radial infinity with the black-hole event horizon. Hence, for multi-transonic accretion, there is no individual existence of physically accept- 
able accretion/wind solution passing through the inner critical (sonic) point, although, depending on the initial boundary conditions, such 
solution may be clubbed with the accretion solution passing through outer critical point, through a standing shock. 

The set [f , A]^ (or more generally [£, A, 7, a]^) thus produces doubly degenerate accretion/wind solutions. Such two fold degeneracy 
may be removed by the entropy considerations since the entropy accretion rates for solutions passing through the inner critical point and 
the outer critical point are generally not equal. For any [f , A, 7, a] £ [f , A, 7, a]^ we find that the entropy accretion rate H evaluated for the 
complete accretion solution passing through the outer critical point is less than that of the rate evaluated for the incomplete accretion/wind 
solution passing through the inner critical point. Since the quantity H is a measure of the specific entropy density of the flow, the solution 
passing through the outer critical point will naturally tend to make a transition to its higher entropy counterpart, i.e. the globally incomplete 
accretion solution passing through the inner critical point. Hence, if there existed a mechanism for the accretion solution passing through the 
outer critical point to increase its entropy accretion rate exactly by an amount 

AH = H(ri")-H(rr), (49) 

there would be a transition to the accretion solution passing through the inner critical point. Such a transition would take place at a radial 
distance somewhere between the radius of the inner sonic point and the radius of the point of inflexion of the homoclinic orbit. In this way 
one would obtain a true multi-transonic accretion solution connecting the infinity and the event horizon, which includes a part of the accretion 
solution passing through the inner critical, and hence the inner sonic point. One finds that for some specific values of [£, A, 7, a]^, a standing 
Rankine-Hugoniot shock may accomplish this task. A supersonic accretion through the outer sonic point (which in obtained by integrating 
the flow starting from the outer critical point) can generate entropy through such a shock formation and can join the flow passing through the 
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inner sonic point (which in obtained by integrating the flow starting from the outer critical point). Below we will carry on a detail discussion 
on such shock formation. 



5.5 The shock formation and related phenomena 

In the present work, the basic equations governing the flow are the energy and baryon number conservation equations which contain no 
dissipative terms and the flow is assumed to be inviscid. Hence, the shock which may be produced in this way can only be of Rankine- 
Hugoniot type which conserves energy. The shock thickness must be relatively small in this case, otherwise non-dissipative flows may 
radiate energy through the upper and the lower boundaries because of the presence of strong temperature gradient in between the inner and 
outer boundaries of the shock thickness. In the presence of a shock the flow may have the following profile. A subsonic flow starting from 
infinity first becomes supersonic after crossing the outer sonic point and somewhere in between the outer sonic point and the inner sonic 
point the shock transition takes place and forces the solution to jump onto the corresponding subsonic branch. The hot and dense post-shock 
subsonic flow produced in this way becomes supersonic again after crossing the inner sonic point and ultimately dives supersonically into 
the black hole. A flow heading towards a neutron star can have the liberty of undergoing another shock transition after it crosses the inner 
sonic point. Or, alternatively, a shocked flow heading towards a neutron star need not encounter the inner sonic point at all. This is because, 
the hard surface boundary condition of a neutron star by no means prevents the flow from hitting the star surface subsonically. 

For the complete general relativistic accretion flow discussed in this article, the energy momentum tensor T^'^, the four- velocity v^, and 
the speed of sound Cs may have discontinuities at a hypersurface E with its normal -q^ . Using the energy momentum conservation and the 
continuity equation, we have: 

[[p«n]'?M=0,[[Tn]'7.=0. (50) 
For a perfect fluid, we thus formulate the relativistic Rankine-Hugoniot conditions as 

[\puTu\] = 0, (51) 
[[TtM^"]] = [[{p + e)vtuV^]]^Q, (52) 
W'^^^n^'n]] = [ [(P + ^)uTl + p] ] = 0, (53) 



where = 1/ \/l — is the Lorentz factor. The first two conditions l l51t and l l52b are trivially satisfied owing to the constancy of the 
specific energy and mass accretion rate. The constancy of mass accretion yields 
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The third Rankine-Hugoniot condition l l53t may now be written as 

7 - 1 \ ^ / cl \ ^ J (7 - c?) + c; 



: 0. (54) 



: 0. (55) 



1 J Vl-l-ci/ I ci(l-«2) 

Simultaneous solution of Eqs. ( I54t and ( I55t yields the 'shock invariant' quantity S( which changes continuously across the shock surface. 
We obtain the analytical expression for 5( in terms of various local accretion variables and initial boundary conditions. The shock location 
in multi-transonic accretion is found in the following way. While performing the numerical integration along the solution passing through 
the outer critical point, we calculate the shock invariant Sh in addition to u, Cs and M. We also calculate Sh while integrating along the 
solution passing through the inner critical point, starting from the inner sonic point up to the point of inflexion of the homoclinic orbit. We 
then determine the radial distance r^h, where the numerical values of Sh, obtained by integrating the two different sectors described above, 
are equal. Generally, for any value of [£, X, 7, a]allowing shock formation, one finds two formal shock locations, one located in between the 
outer and the middle sonic point and the other one located in between the inner and the middle sonic point. The shock strength is different for 



the inner and for the outer shock. According to the standard local stability analysis l lYang & KafatosI ( Il995h ). for a multi-transonic accretion, 
one can show that only the shock formed between the middle and the outer sonic point is stable. Hereafter, whenever we mention the shock 
location, we always refer to the stable shock location only. 

We find that the shock location correlates with A. This is obvious because the higher the flow angular momentum, the greater the 
rotational energy content of the flow. As a consequence, the strength of the centrifugal barrier which is responsible to break the incoming 
flow by forming a shock will be higher and the location of such a barrier will be farther away from the event horizon. However, the shock 
location anti-correlates with £ and 7. This means that for the same £ and A, in the purely non-relativistic flow the shock will form closer 
to the black hole compared with the ultra-relativistic flow. Besides, we find that the shock strength anti-correlates with the shock location 
Tsh, which indicates that the closer to the black hole the shock forms , the higher the strength and the entropy enhancement ratio are. The 
ultra-relativistic flows are supposed to produce the strongest shocks. The reason behind this is also easy to understand. The closer to the 
black hole the shock forms, the higher the available gravitational potential energy must be released, and the radial advective velocity required 
to have a more vigorous shock jump will be larger. Besides we note that as the flow gradually approaches its purely non-relativistic limit, 
the shock may form for lower and lower angular momentum, which indicates that for purely non-relativistic accretion, the shock formation 
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may take place even for a quasi-spherical flow. However, it is important to mention that a shock formation will be allowed not for every 
[£, A, 7, a] G [£, A, 7, a]^, The numerical value for <S( will be the same along both the accretion through the outer and the inner sonic point, 
only at certain points (the shock locations) for a specific subset of [£, A, 7, a]^, for which a steady, standing shock solution will be found. 

We further find that the shock location correlates with the black hole spin parameter, whereas the shock strength and the entropy 
enhancement ratio at the shock anti correlates with the Kerr parameter. The post to pre shock velocity correlates with the black hole spin, 
whereas the post to pre shock temperature, density and pressure anti correlates with the Kerr parameter. 



5.6 The 'asymmetrically hunched fish' and the heteroclinic orbits 

As we further increase A, the shock location recedes outward, and the shock strength as well as Rj. decreases monotonically. Meanwhile, the 
homoclinic path embracing the middle critical point starts filling up the entire region accessible to it. Finally, when A approaches the value 
A5~3. 57401752, the interface between the A and the Ai regions is approached, and one obtains the characteristic phase portrait as that of 
panel (v). 

The above mentioned line of interface, as has already been discussed, represents the state where the entropy accretion rate for solution 
passing through the inner critical point is exactly the same for the solution passing through the outer critical point. Integral curves on the phase 
portrait are allowed to intersect only at the critical points. Hence two saddle points in this case will be connected through heteroclinic orbits. A 
heteroclinic orbit is generally a trajectory on a two dimensional phase portrait which connects one saddle type equilibrium points to the other 
saddle type equilibrium point. If the two saddle points are the same, then it becomes a homoclinic orbit. Existence of the heteroclinic orbit 
is associated with the stability issues in the dynamical systems of a two dimensional vector field (velocity filed, for example). A heteroclinic 
orbit is contained in the stable manifold of one saddle point and in the unstable manifold of the other saddle point. A fundamental fact arising 
from the theory of the structural stability is that if two saddle type critical points are connected through the heteroclinic orbit, then the local 
phase portrait near this orbit can be changed by an arbitrarily small smooth perturbation. In effect, a perturbation can be chosen such that, 
in the phase portrait of the perturbed vector field, the saddle connection is broken. Thus, in particular, a vect or field with tw o saddle points 
connected by a heteroclinic orbit is not structurally stable with respect to the class of all smooth vector fields 1 Chicone ( 20061) ). 



Drawing the phase portrait for the above mentioned topology will entail the tuning of the flow parameters (and the initial boundary 
conditions) with infinite numerical precession. The solution topology looks like a fish, which has an asymmetrically placed hunch on its 
back. The collection of all the [£ — A] pair, for which such 'fish topologies' are produced, is thus a line characterizing the condition of the 
heteroclinicity. 



5.7 Topology from region Ai and the pair annihilation of the critical points 

A further increase of A allows to enter in the Ai region. A representative diagram for such flow topology is shown in the panel (vi). The 
converse behaviour (compared to the multi-transonic accretion) of the phase portrait for multi-transonic accretion) is quite obvious, here the 
governing principle, as already been mentioned in Sect. 14.31 is that the entropy accretion rate pertaining to the inner critical point is greater 
than that of the outer critical point, and hence the homoclinic orbit passes through the outer critical point. 

The topological flow profile for the present case is different from that of the multi-transonic accretion. Here the homoclinic orbit passes 
through the outer critical point. The procedure for finding the complete flow topology with shock is essentially same as that for finding the 
complete topology of a true multi-transonic shocked accretion flow as mentioned earlier. To find the shock location, here one calculate the 
numerical values of Sh along the wind solution passing through the inner critical point are compare with the numerical values of Sh along 
the wind solution passing through the outer critical point, and the stable shock location is found accordingly. Extremely strong shocks are 
found to form very close to the black hole (on the off equatorial plane) for this situation. We will show in our future work that such wind 
solutions with strong shocks are useful in explaining the GRB light curves (Janiuk, Czemy & Das, in preparation). 

For higher value of A, the homoclinic orbit through the outer critical point start shrinking, and finally disappear with pair-annihilation 
of the saddle type outer critical point and the centre type inner critical point once we cross the hypersurface separating the Ai and the mono- 
transonic accretion (passing only through the inner critical point) zone. All that is left behind is a single saddle type critical point, through 
which mono-transonic accretion and wind passes. 



6 CONTINUITY IN THE PARAMETER SPACE AND THE HYSTERESIS EFFECT 

The heteroclinic fish like solution as described in the section 5.6 marks the transition between two different flow topologies, a multi-transonic 
accretion with a built-in homoclinic orbit (as shown in the successive panel figures (ii) - (iv)) and a mono-transonic accretion with an 
additional non accessible homoclinic orbit (as shown in the panel figure (iv)). It is to be understood that the solutions with shock for the value 
of the flow angular momentum higher than A5 are of completely different character compared to the solution parametrized by the value of 
angular momentum lower than A5. We illustrate this issue through the following procedure: 

We set rio = lOr-g as our reference point. For a transonic solution, we then calculate the density of the accerting matter at rio as a 
function of A, keeping the other initial boundary conditions [£,7, a] fixed. This has to be done in several steps. For mono-transonic accretion 
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Figure 3. In the figure, Prout Pi-i,, represent the density of the accreting material measured at lOr'g along the transonic accretion passing through the 
outer and the inner sonic points, respectively. The long dashed line coiresponds to the range of the angular momentum for which a steady standing Rankine 
Hugoniot shock forms. The figure has been drawn for [£ = 1.0000033, 7 = 4/3, a = 0.3]. 



passing through the outer sonic point only, we can have only one choice to calculate p at rio, that is to calculate the density at rio only on 
the solution passing through the outer sonic point, Tout, and denoting it by praut- We do it without any ambiguity for a span of value of 
A, starting from a very small value (almost equal to zero, just when the departure from the pure Bondi type accretion stars) upto the value 
A = Ai. Beyond Ai, for A > Ai, in principle there are two choices for the density calculation at rio. One is along the flow passing through 
rout, to obtain praut for the multi-transonic flow. Another option is to compute the density pri„ along the multi-transonic accretion flow 
passing through the inner sonic point ri„. Both pra^t ™d p^^ can be calculated as a function of A for the span Ai < A < A5. However, 
since for the span Ai < A < A3 shock does not form, pr^^ is a meaningless quantity to define for this span of A, and thus for any value of the 
angular momentum upto A = A3, density at rio means pr^ut only. However, the situation is different for A3 < A^ A5, where both Pr^ut and 

can have feasible meaning. However, praut is practically meaningful if shock does not form (since it has been thoroughly verified that 
the lowest possible shock location j-^'""^"'" for the entire range of our choice of [£, A, 7, a] is always greater than rio), and pr^^ is practically 
meaningful only if the shock forms. Beyond A5, for A > A5, once again one has mono-transonic accretion (with an additional homoclinic 
orbit of course), now passing through the inner sonic point Tin. Hence the measured density at rio for A > As implies pr^^ only. 

Hence, we have p,.„„t defined over the range for any value of the flow angular momentum upto A = A5, and, on the contrary, pr^^ is 
defined only over the range A^A3. In figure 3, we plot, in logarithmic scale, p,-„„t using the solid line (red coloured in the online edition) 
and pr;,, using the solid line (in the region IV) and the dashed line (in the region III) respectively, as a function of the specific flow angular 
momentum A. 

The local density changes by a factor over 500 at A5 which well illustrates the global quantitative change in the flow. However, when we 
add the solutions passing through the inner and the outer sonic points with a shock to the same plot, those solutions form a continuation of the 
previous branch below A5. The shock for A just below A5 is extremely weak (with the corresponding shock strength having value very close 
to the unity), and located infinitesimally close inward of the outer critical point. With a further decrease of A the shock strength measured as 
a ratio of the Mach numbers above and below the shock increases monotonically, and the shock position moves inward. 

The branch with a shock also does not cover the whole range of angular momentum values - there is a minimum value of A (for a given 
Bernoulli constant, the polytropic index and the black hole spin) for which the Raiikine-Hugoniot conditions (see Sect. 15. 5t can be satisfied. 
This value, marked as A4 in the section IsTtI describes the discontinuous change in the solution along the branch with a shock with a further 
drop in the angular momentum, as for still lower values of A only the solutions without shock exist. The local density change in the flow at 
rio accompanying the slight change of angular momentum along this branch is by a factor 4, much less than the factor 500 at the transition 
at As for the shock-free solutions. 

Since the system is likely to choose a solution close to the previous one under the slight change of the outer boundary conditions (here, 
the angular momentum), the choice between the solution with a shock and without a shock for A4 < A < A5 will depend on the previous 
history of the system. If the initial value of the angular momentum is lower than A4 and slowly increases the system will follow the lower 
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Figure 4. The variation of the shock strength (sohd hne, red coloured in the online version) and the entropy accration rate ratio at the shock (dotted line, blue 
coloured in the online version) as a function of the flow angulai' momentum for a fixed set of initial boundary conditions [E = 1.0000033, 7 = 4/3, a = 0.3] 
Shock strength anti-correlates with the shock location. As a result, shocks formed closer to the black hole are stronger 



branch (p, „„t) in Fig. 3 without a shock and will show a dramatic change when passing from A < A5 to A > A5. If, on the other hand, the 
system evolution starts with A > A5 and the angular momentum of the infalling material slowly decreases the system is more likely to follow 
the upper branch in Fig. 3, without any discontinuity in the appearance down to A4 but the solution will display an increasingly stronger 
shock with a quasi-stationary increase of the angular momentum of the matter. 

We thus expect a hysteresis effect in the system behavior under the conditions of slowly varying initial angular momentum of the 
accreting matter. This may in turn lead to some asymmetry in the lightcurves of the systems accreting in such a mode. The local density 
at lOr-g is an indicator of the system luminosity (the synchrotron emission increases with the plasma density and magnetic field). The 
increase in the local density indicates an increase in a source brightness. Since the discontinuity in the density profile is larger when the 
angular momentum and consequently the density increases (at A5) the brightening should be faster/stronger than the subsequent faintenning 
of the source (evolution along the upper branch (pri„) in Fig. 3, down to A4). However, detailed study of this effect and comparison with 
observational data is beyond the scope of the present paper. 

In this context, it is to be noted here that the shock strength, which correlates with the entropy enhancement ratio at the shock, 
monotonically increases with decreasing A, and at the same time anti-correlates with the shock location itself. The variation of the shock 
strength (solid line marked by M_ /M^) and (dotted line marked by R^) with A is shown in the figure 4. The strongest shock thus forms 
closer to the black hole. From observational data, it is expected that the QPO profile of the SgrA* black hole corresponding to a small length 
scale shock formation. Shocks formed for the angular momentum A~A4 thus, might corresponds to such QPO behaviour. 



7 A LYAPUNOV LIKE TREATMENT 

As we already have observed, for an extended region of the parameter space describing the multi-transonic accretion flow, shock does not 
form in spite of the the availability of two saddle type sonic points to the accretion flow to pass through. This region has been marked as Ans, 
and is bounded by the span of the specific angular momentum [A2 — A3] . One might raise the question about what dictates the flow not to 
make a shock transition even if the shocked flow can have the saddle type inner sonic point available to pass through. In other words, given all 
the favourable conditions, i.e., the presence of two real physical sonic points, and a larger than one entropy accretion rate ratio for the flows 
passing through these sonic points (with higher entropy accretion rate for the flow passing thorough the outer sonic point), what prohibits 
the Rankine Hugoniot condition to be satisfied, as it gets satisfied for the region As? In reality, no information can readily be obtained from 
the initial boundary conditions [£, A, 7, a] as well as the measure of the entropy accretion rate (which, of course, is purely a function of 
[f , A, 7, a], and hence, does not impose any additional condition) which may provide some certain well defined conditions predicting about 
the validity of the Rankine Hugoniot conditions over certain region of the parameter space describing the multi-transonic accretion flow. 
Usual trends in the contemporary literature 



Chakrabartilll989l : Da3l2002 ) and references therein) to explain the origin and the charac- 
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teristic features of the Ans region is to make hypothesis that the shock may become oscillatory in that region, in absence of the fine tuning 
condition of the entropy generation at the shock (that the difference between the entropy accretion rate should exactly be equal to the excess 
amount of H produced at the shock), and the region Ans thus has been attributed to ex plain, in analogy with the results obtained from the 
full numerical time dependent simulation of the shocked axisymmetric flow (see, e.g., i Sponholz & Molteni 1994; Okuda et al J 2007 1 and 
references therein), the formation of the quasi-periodic oscillation of the black hole candidates. This approach, however, is fundamentally 
inconsistent from the mathematical perspective. The parameter space corresponding to the formation of the Ans region (of any region of 
the parameter space shown in the figure 1) is obtained from a completely stationary configuration, where the mother equations have been 
constructed as a set of explicit time independent equations. The notion of the oscillation of the shock (or the very phrase ' 'oscillation' itself) 
is associated with the dynamical temporal evolution of the configuration governed by a set of full time dependent equations. Hence, although 
for full time dependent numerical simulations, the appearance of the shock oscillation can be considered as a consistent physical scenario, 
such shock oscillation prescriptions can never be applied to clarify any results which has been obtained by solving a set of explicit steady 
state equations. 

The above issue prompted us to put forward a theoretical condition, which, within the framework of the stationary configuration, will 
be able to explain why the entire region of the multi-transonic accretion parameter space will not allow a stationary shock solution, and, 
instead, only a restricted part of the region will produce a shocked accretion flow. We propose a Lyapunov like treatment of the transonic 
accretion solution, within a complete stationary framework, to accomplish the above task. It is not quite uncommon in the literature to 
reply on a Lyapunov like treatment, within the configuration of a stationary set up which contains an equilibrium point (the critical point) 
to study the state transition of such systems, subjected to a small variation of the overall parameter space. Such a methodology have been 
applied to several systems, f or example, to describ e the state transition of a spin galas phase in terms of the T = Q critical point where T 
is the absolute temperature ( dBrav & Moordl 19871) ). or to study the thermod ynamic stability of the p rotein like heteropolymers against the 
biological mutation process, by considering such process as a perturbation ( dBussemaker et al.lll997n ). In what follows, we propose a one 
dimensional effective Lyapunov index Leff. Leff will exclusively predict which final state - a shocked flow passing through both the outer and 
the inner sonic points, or a shock free flow passing only through the outer sonic point - a multi-transonic accretion will follow, subjected to 
the specific initial boundary conditions, and the infinitesimal variation of the control parameter describing the secular evolution of the flow. 

Proposition 7.1. On a two dimensional parameter space spanned by [£, A] as shown in the figure 1. let us consider two transonic solutions 
characterized by infinitesimally different values of the specific flow angular momentum \i and A/, satisfying the condition 

Lt.^o |A, -Ay l=6 (56) 

Let u\ (r) and U2 (r) be the values of the advective velocity at a certain value of the radial distance r for two transonic solutions characterized 
by very close values of the angular momentum Ai and A2 respectively, with £ and 7 kept fixed. Let u\{r — R) and U2{r — R) be the 
corresponding velocities for the above-mentioned transonic flow solutions at the radial distance r — R, with the obvious requirement that 
R < r. We then define the effective Lyapunov exponent Leff, so that the following condition holds: 

uijr - R) - U2{r - R) ^ ^L„ffl ^^^^^ 
ui{r) ~ U2{r) 

The effective index Leff can also be calculated by varying £ instead of A. 

Without any loss of generality, one can normalize the constant to appear in front of the e'"'"^^ and can define the effective Lyapunov 
exponent as: 

^^^^^ 1 r^L-^l^n,^! (58) 



R { ui(r) — U2{r) 

It is obvious from the expression of Leff that for accretion, the effective Lyapunov index is a measure of how drastically the output (Mach 
number variation in the flow topology) will change with respect to an arbitrary small perturbation in the control parameter (the specific 
angular momentum A, for our case). Hence there is no harm in assuming that Leff does share the properties of a Lyapunov exponent in the 
context tha t a higher value of the m odulus of Leff indicates a greater available volume of the phase space (in connection to the Kaplan- Yorke 
conjecture ( Kaplan & YorkellT979h ). 

We now consider two different values of angular momentum infinitesimally separated and lying at two different side of the value A = Ai 
(the value of the angular momentum corresponding to the boundary of the O and the A region of the parameter space diagram as shown in 
the figure 1), as A± where A_ = Ai — e and A+ = Ai + e for e^O. For A = A_, one has a mono-transonic accretion only through the outer 
sonic point, and A = A+ provides the multi-critical accretion. We now calculate the Leff for [r — lOrj, R = 0.2rj, A_ = 2.839, A+ — 2.84], 
where u\ corresponds to the solution characterized by A_ and U2 corresponds to the solution characterized by A = A2, all for the initial 
boundary condition [£ — 1.0000033, 7 = |, a = 0.3] for the following two different situations as described below: 

(i) For A — A+, we fist assume that the accretion flow, even formally being multi-transonic, passes only through the outer sonic point 
(which is the reality because solution characterized by A+ is a shock free solution from Ans). We then calculate U2{r) and U2{r~R) for flow 

fxout \outl 

passing through the outer sonic point. For such situation, we denote the effective Lyapunov index as Leffl^ ^ ' + J, where it is understood 
that the mono-transonic accretion characterized by A_ passes only through the outer sonic point and the the effective Lyapunov index is 
being calculated at r = lOrg for accretion solution characterized by A+, which has been considered to pass through the outer sonic point. 
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where the outer sonic point is larger than r = lOrg. Within the bracket. A?"* has been written first, and then . This sequence indicates 
that the change of the specific angular momentum has been made such a was that a transition is being induced from the mono-transonic 
accretion (passing through the outer sonic point) characterized by A_ to the multi-transonic accretion (but here again, passing only through 
the corresponding outer sonic point, and does not encounter the inner sonic point because shock does not form in this case) characterized by 
A+. The modulus of the numerical value of the effective Lyapunov index, which has been denoted by mod ^Lefff'*'" ''^+ 1^, is found to be 
~ 1.6631 for the present situation. 

(ii) Here for A+, we assume, for a hypothetical situation, what would have been the value of the effective Lyapunov index at r = lOrg, 
if, for whatever reason, the muli-critical accretion would finally pass through the inner sonic point. Hence U2 (r) and U2 (r — R) have been 
calculated for flow passing through the inner sonic point, and the corresponding effective Lyapunov index LefJ"^- '^+1 has been computed. 
The the corresponding modulus of the numerical value of Leffl^^ ''^+ \ which is denoted as mod ^Leffl^^ '^+1^ is found to be ~1. 60738. 

We thus clearly have 

mod (^Lefft^-"''^""]) > mod (lefft^-" *'^+]) (59) 

We now apprehend what does the above inequality indicate. 

If one makes a transition from the mono-transonic flow to the multi-transonic flow by an incremental smooth change of the specific 
angular momentum, Leff is a essentially a measure of the probability of the outcome of final state corresponding to such a transition. The 
effective Lyapunov exponent has been defined in such a way that the above findings allows us to make the following statement: 

Conjecture 7.2. For a certain transonic accretion state characterized by a control parameter d (Ci = A_ = Ai — e for the present purpose), 
if there are two possible transitions to two different transonic states, both characterized by the same value of the final control parameter Cf 
(cf = A+ = Ai + e for the present purpose), differing from Ci by an infinitesimal measure e— ^-O, and if such transitions can be quantified 
by two different effective Lyapunov indices Ljff and L^ff, then the transition characterized by the higher absolute value (the modulus of the 
numerical value) of the effective Lyapunov index is a more probable outcome with greater likelihood, since a higher value of the Lyapunov 
index indicates a greater available occupancy of the phase volume for the solution. 

This means that according to the inequality l l59t . the probability of finding (at r = lOr^) the accretion flow passing through the outer 
sonic point (shock free multi-criticality) is higher than the probability of finding it to pass through the inner sonic point (shocked multi- 
transonic accretion). And that is exactly what is observed for the Ans region. We calculate, for a wide range of r and R, with r lying close 
to the inner sonic point, that is, closer to the black hole in comparison to any 'anticipated' artificial shock location, so as to calculate the Leff 
for an 'artificial' post shock flow (had it been the case that the shock would, at all, form, which clearly does not, in reality). This procedure 
guarantees the consistency of checking the probability of two possible outcomes at a length scale where one can distinguish a shocked and a 
shock free solution. For all possible configurations, we consistently found that the inequality l l59t always holds good. 

What about the true multi-transonic accretion flow structure allowing a standing shock? We can show that sufficient room for such 
flow configuration is available as well. For this purpose, consider two infinitesimally separated value of A, A± at the opposite sides of the 
boundary separating the A and the Ai region. Here A_ — X5 ~ e provides a multi-transonic accretion and A-|- = A5 + e provides a transonic 
accretion passing only through the inner sonic point (and associated with a homoclinic orbit characterized by an additional pair of critical 
point, none of which a real physical accretion solution can pas through). We now calculate the corresponding effective Lyapunov indices for 
a decremental smooth change of the angular momentum. The unique initial state (accretion passing only through the inner sonic point) is 
characterized by A-|-. Two final states, both characterized by A_, are available. The first option for such final state is a hypothetical shock 
free multi-critical accretion passing through the outer sonic point only, on which U2{r) and U2{r — R) are measured at some specified value 

r\in \Outl 

of r and R. The corresponding effective Lyapunov index being Leff ^ + ' - J . For the other variant (the actual case), the shock forms and the 
accretion finally passes through the inner sonic point. The corresponding index is LefF[^+ '^-^ 

A suitable choice of r is trivial, it has to be smaller than the smallest possible shock location for A_ type flow, so as to the length scale 
for our interest lies in the post-shock region. All possible variation of such r and R confirms the generic nature of the following inequality 

mod (^Leff[^+'^-l) > mod (^L - eff[^+'^+*]) (60) 

Lemma 7.3. For an incremental smooth change of the flow specific angular momentum, a transition from a mono-critical mono-transonic 
flow to a multi-critical flow will not produce a shock, however, for a decremental smooth change of A, a transition from a mono-transonic 
accretion to a multi-critical accretion must accompany a standing shock, and the multi-critical accretion is a true multi-transonic one. 

It is interesting to note that the above lemma has one to one correspondence with our findings related to the hysteresis effect as described 
in the previous section. 

We thus define an independent quantity, within a framework of a purely stationary configuration, namely the effective Lyapunov index 
Leff, which, for the first time in literature, can analytically explain why one obtains shock free solutions even for a multi-critical flow 
configuration. 
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8 DISCUSSION 

In the present paper we study the topology of the low angular momentum flow onto a black hole as a function of the flow global parameters. 
Such a study may indicate the behaviour of the time evolution of the flow seen as a sequence of quasi-stationary stages. 

For a given set of other parameters (the Bernoulli constant, adiabatic index and the Kerr parameter) the change in the angular momentum 
of the flow (determined by the outer boundary conditions) leads to a sequence of solutions passing through two discontinuous changes in the 
flow topology, and consequently, in the flow local parameters. 

The discontinuity at A5 corresponds to a heteroclinic solution (a 'fish' topology) with a flow line connecting the inner and the outer 
critical points. As we argue in Sect.|6l the physical flow is likely to change discontinuously at A5 if the angular momentum increases from 
A < A5 to A > A5 as a result of systematic secular evolution, and then solutions both for A < A5 and A < A5 do not develop a shock. 
However, if in the course of secular evolution the angular momentum decreases from A > As to A < As the systems evolves through As 
without showing any discontinuous change but instead develops a shock, with the strength of the shock increasing with the As /A ratio. 

The above mentioned finding has further been supported by a Lyapunov like treatment of the state transition of the flow, as discussed in 
detail in Sect.|7] 

The discontinuity at A4 corresponds to the minimum value of the angular momentum for which the solutions with shocks exist. It is 
likely to show up if the angular momentum systematically decreases, first from A > As to A < As (following the branch with shocks) and 
then down below A < A4. The solution must then change discontinuously to the one without a shock. The reverse evolution does not lead to 
discontinuity as the flow then proceeds as a shock-free solution continuously up to As where it reaches the critical 'fish' topology. 

The analytical study of the solution stability close to 'fish ' topology (A oc As) cannot be performed since infinitesimal changes in the 
solution parameters lead to the arbitrarily large change in the flow parameters. The problem should be addressed through time-dependent 
dynamical simulations. 

The issue of the topological change in the flow is generic and not related to any specific choice of the other flow parameters. We will 
show in our next work (Das & Czerny, in preparation) that the 'fish' topology always appears as a boundary between the inflow through 
the inner critical point and through the outer critical point, irrespective of the flow geometry (as long as it contains some non zero angular 
momentum), the specific metric, and the choice of the equation of state. 

Our discussion shows that both solutions with, and without a shock, are likely to be met in real sources accreting material with low 
angular momentum, and the actual choice of the solution by the system depends on the flow history. 

If the accretion parameters in an astrophysical system are in the range somewhat broader than A4 - As and vary in time the system 
will show sudden rapid changes in the luminosity accompanying the sudden flow restucturization. We qualitatively predict the appearance 
of flares with more rapid increase of the luminosity (when crossing from A < As to A > As) followed by a slow decay (evolution partially 
along the shock branch down to A < A4). Quantitative estimates would require computations of the flow emissivity. 

The appearance of the homoclinic orbit on the phase portrait is not unique for the axisymmetri c rotating flows on l y. For spherically sym- 



metric accretion with self gravity and radiative pr ocesses incorpora ted into it Karkowski. Malec & Roszkowski 1 20081) : Karkowski. Malec, Roszkowski & Swie 



1 20091) . such homoclinic paths would also appear] iRembiaszl |2009|) . making the spherically symmetric accretion to be a bi-critical flow. It ap 



pears that the emergence of a homoclinic orbit is due to the smooth decrement of the radial inflow velocity, due to the centrifugally supported 
barrier for rotating flow, while due to the radiative losses for the spherical accretion. 

In this work, we have used a relatively lower value of the Kerr parameter (a = 0.3) as a representative one. This does not, by any means, 
limit the applicability of our methodology for the case of black holes possessing higher value of the spin angular momentum because the 
Kerr parameter has been treated as a free parameter in our construction and any value of a can easily be plugged in into our model set up. For 
example, we obtain consistent shock location for the value of a as high as 0.99 (ad for values of 7 other than 4/3) as well. We can, compute 
the shock location and all other shock related quantities in general as a function of the Kerr parameter as well, and that has been presented in 
some other work (Barai et al. in preparation). 

Some of the results presented in this paper may also apply to ORB. In this event, modeled as a failed supernova, the star first collapses due 
to the fuel exhaustion, the collapse is halted close to the core (at ~ 500km) and replaced by expansio n. Part of the stellar envelope is rejected 



but the neutrino cooling stops part of the material and subsequently the slow fallback stage takes place lMacFadven. Wooslev & Hegeg d2001l) 
onto a newly formed compact object. This accretion stage of the long gamma-ray burst can be described within the frame of our model 
assuming that the angular momentum of the back falling envelope is not large. If the accretion flow proceeds through a solution with a 
shock, the sho ck-related pe r iodic oscillations may in a natural way explain for example the fast optical variability of the Naked-Eye burst 
GRB080319BlBeskin et all J2OO9I). 
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